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Q . Abstract 

<!-> ■ 

The q-state Potts model can be defined on an arbitrary finite graph, 

and its partition function encodes much important information about 
that graph, including its chromatic polynomial, flow polynomial and 
reliability polynomial. The complex zeros of the Potts partition func- 
tion are of interest both to statistical mechanicians and to combina- 
C^ , torists. I give a pedagogical introduction to all these problems, and 

^ ' then sketch two recent results: (a) Construction of a countable family 

I , of planar graphs whose chromatic zeros are dense in the whole complex 

' ^ ' g-plane except possibly for the disc |g — 1| < 1. (b) Proof of a universal 

upper bound on the g-plane zeros of the chromatic polynomial (or an- 
fj ' tiferromagnetic Potts- model partition function) in terms of the graph's 

maximum degree. 



(^ '. 1 Introduction 

O. 

*/^ ■ The Potts model [||, y, ^ plays an important role in the general theory of 

critical phenomena, especially in two dimensions Q, ^, |^, and has applica- 

Q^ ' tions to various condensed-matter systems Q- Ferromagnetic Potts models 

^^ ■ have been extensively studied over the last two decades, and much is known 

about their phase diagrams |2|, ^ and critical exponents fS], ^ |^. But for 

antiferromagnetic Potts models, many basic questions remain open: Is there 

a phase transition at finite temperature, and if so, of what order? What is 

the nature of the low-temperature phase(s)? If there is a critical point, what 

Q I are the critical exponents and the universality classes? The answers to these 

J-J ' questions are expected to be highly lattice-dependent, in sharp contrast to 

the universality typically enjoyed by ferromagnets. This suggests studying 

k^ I antiferromagnetic Potts models on general graphs, in order to learn which 

^ ' combinatorial properties of the lattice correlate with which properties of the 

phase diagram. 
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So let G = (y, E) be a finite undirected graph with vertex set V and edge 
set E. For each positive integer q, let Pg{q) be the number of ways that 
the vertices of G can be assigned "colors" from the set {1, 2, . . . ,q} in such 
a way that adjacent vertices always receive different colors. It is not hard to 
show (see below) that Pg{q) is the restriction to Z^ of a polynomial in q. 
This (obviously unique) polynomial is called the chromatic polynomial of 
G, and can be taken as the definition of Pg{q) for arbitrary real or complex 
values of q^ 

The chromatic polynomial was introduced in 1912 by Birkhoff [ll[|. The 
original hope was that study of the real or complex zeros of Pg{q) might 
lead to an analytic proof of the Four-Color Conjecture, which states that 
Pg(4) > for all loopless planar graphs G. To date this hope has not been 
realized, although combinatoric proofs of the Four-Color Theorem have been 
found ||T2| , p!3[| . Even so, the zeros of Pg{q) Biie interesting in their own 
right and have been extensively studied. Most of the available theorems 
concern real zeros, but there has been some study (mostly numerical) of 
complex zeros as well, by both mathematicians and physicists (see [|14| for 
an extensive list of references) . 

A more general polynomial can be obtained as follows: Assign to each 
edge e G £" a real or complex weight Ve- Then define 



ZGiqdve}) = E n 



'^ + VeS{(7xiie),(^X2ie)] 



fl.l) 



where the sum runs over all maps a: V ^ {1, 2, . . . , q}, the 6 is the Kronecker 
delta, and xi{e),X2{e) E V are the two endpoints of the edge e. It is not 
hard to show (see below) that ZG{q,{ve}) is the restriction to q G Z+ of 
a polynomial in q and {ve}- If we take Ve = —1 for all e, this reduces to 
the chromatic polynomial. If we take Ve = v for all e, this defines a two- 
variable polynomial ZG{q, v) that was introduced implicitly by Whitney |p!5|, 



16, 17 1 and explicitly by Tutte [|^, 19 1; it is known variously (modulo trivial 



changes of variable) as the dichromatic polynomial, the dichromate. 



the Whitney rank function or the Tutte polynomial |gO|, |21[| . 



In statistical mechanics, (1.1) is the partition function of the q-state 



Potts model ||T|, |2|, ^ with Hamiltonian 



eeE 



^ See M, H] for excellent reviews on chromatic polynomials, and [tldl for an extensive 
bibliography. 



here we have written 

Ve = e^-^-^ - 1 . (1.3) 

A couphng Je (or Ve) is called ferromagnetic if Jg > {ve > 0) and 
antiferromagnetic if — oo < Je < (— 1 < fg < 0). The chromatic polyno- 
mial thus corresponds to the zero-temperature limit of the antiferromagnetic 
Potts model. 

To see that ZQ{q,{ve}) is indeed a polynomial in its arguments (with 
coefficients that are in fact or 1), we proceed as follows: In (|l.lj), expand 
out the product over e a E, and let E' C E he the set of edges for which 
the term Vp6a , , „ , a is taken. Now perform the sum over configurations 
{(Tx}'- in each connected component of the subgraph (y,E') the spin value 
ax must be constant, and there are no other constraints. Therefore, 

ZGiq,{ve}) = Yl l'^""'^ n^-' (1-4) 

E'CE ee-B' 

where k{E') is the number of connected components (including isolated 
vertices) in the subgraph {V,E'). The expansion ( |1.4p was discovered by 
Birkhoff [^] and Whitney |jl5| for the special case Ve = —1 (see also Tutte 
|l8| , p^); in its general form it is due to Fortuin and Kasteleyn |^, ^ 
(see also [^). We take (|l.4|) as the definition of Zciq, {ve}) for arbitrary 
complex q and {ve}- 

Let us note that a vast amount of combinatorial information about the 
graph G is encoded in Zdq, {^e})- Special cases of Zq include not only the 
chromatic polynomial (v = —1) but also the fiow polynomial (v = —q), the 



reliability polynomial (q = 0) and several other quantities [20|. 

In statistical mechanics, a very important role is played by the complex 
zeros of the partition function [^. Recall that a phase transition occurs 
whenever one or more physical quantities (e.g. the energy or the magneti- 
zation) depend nonanalytically on one or more control parameters (e.g. the 
temperature or the magnetic field). Now, such nonanalyticity is manifestly 



impossible in (l.l)/( |l.4D for any finite graph G. Rather, phase transitions 
arise only in the infinite- volume limit. That is, we consider some countably 
infinite graph Goo = (V'oc-E'oo) — usually a regular lattice, such as U^ with 
nearest-neighbor edges — and an increasing sequence of finite subgraphs 
Gn = iyn^En)- It cau then be shown (under modest hypotheses on the Gn) 
that the (limiting) free energy per unit volume 

IgM^v) = hm \Vn\-^\og ZgM.v) (1.5) 

n— >oo 

exists for all nondegenerate physical values of the parameters, namely either 



(a) q integer > 1 and — 1 < u < oo [using (|1.1|): see e.g. [26, 
Section 1.2]] 
or (b) q real > and < t) < cxo [using ([L.4|): see [^, Theorem 
4.1] and M, M]. 



This hmit /^^ iQ,v) is in general a continuous function of v; but it can fail to 
be a real-analytic function of v, because complex singularities of log Zg„ {q, v) 
— namely, complex zeros of ZG„{q,v) — can approach the real axis in the 
limit n ^ oo. Therefore, the possible points of physical phase transitions 
are precisely the real limit points of such complex zeros |2^. As a result, 
theorems that constrain the possible location of complex zeros of the par- 
tition function are of great interest. In particular, theorems guaranteeing 
that a certain complex domain is free of zeros are often known as Lee- Yang 
theorems. 

2 A Lee- Yang Theorem for Chromatic Polynomi- 
als? 

Let me now review some known facts about the real zeros of the chro- 
matic polynomial Pciq), in order to motivate some conjectures concerning 
the complex zeros: 

1) It is not hard to show that for any loopless graph G with n vertices, 
{—!)"' Pciq) > for real 9 < [§. It is then natural to ask whether the 
absence of negative real zeros might be the tip of the iceberg of a Lee- Yang 
theorem: that is, might there exist a complex domain D containing (— oo, 0) 
that is zero-free for all Pq? One's first guess is that the half-plane Reg < 
might be zero-free [|3^. This turns out to be false: examples are known of 
loopless graphs G whose chromatic polynomials have zeros with real part as 



negative as ss —0.7 [31, 32, 33, ^, |3^, |3^, 37[. Nevertheless, it is not ruled 
out that some smaller domain D D (— cxd,0) might be zero-free. 

2) For any loopless planar graph G, Birkhoff and Lewis p^ p roved in 
1946 that Pclq) > for real g > 5; we now know that Pg(4) > Op], [l|[; and 
it is very likely (though not yet proven as far as I know) that Pg{q) > also 
for 4 < q < 5. Thus, it is natural to conjecture that might exist a complex 
domain D containing (4, oo) [or (5, oo)] that is zero-free for all planar Pq- 
One's first guess might be that Reg > 4 works. This again turns out to 
be false: examples are known of loopless planar graphs G whose chromatic 
polynomials have zeros with real part as large as « 4.2 [^, |3^. Nevertheless, 
it is not ruled out that some smaller domain D D (4, oo) might be zero-free. 

As with most of my conjectures, these two are false; but what is inter- 
esting is that they are utterly, spectacularly false, for I can prove: 



Theorem 2.1 |39| There is a countahly infinite family of planar graphs 



whose chromatic zeros are, taken together, dense in the entire complex q- 
plane with the possible exception of the disc \q — 1\ < 1. 

The graphs in question are "theta graphs" @s,p obtained by parahel-connecting 



p chains each of which has s edges in series. Theorem 2.1 is in fact a corollary 



of a more general result for the two-variable polynomials Zciq^v): 



Theorem 2.2 |39| Fix complex numbers qo,Vo satisfying \vq\ < \qQ +fo|. 
Then, for each e > 0, there exist numbers sq < oo and po{s) < oo such that 
for all s > sq and p > po{s): 

(a) Zq^^{ ■ ,vo) has a zero in the disc \q — qo\ < e. 

(b) ZQ^^{qQ, ■ ) has a zero in the disc |t; — fo| < e. 



(Setting vq = —1, Theorem |2.2| (a) implies Theorem 

The intuition behind Theorem ^^ is based on recalling the rules for 
parallel and series combination of Potts edges: 

Parallel: Vejj =1^1 + ^2 + ■yi'y2 (mnemonic: 1 + v multiplies) 

Series: Vejj= viV2/{q + vi+V2) (mnemonic: t7/(g + w) multiplies) 

In particular, if < \v/{q + v)\ < 1, then putting a large number s of edges 
in series drives the effective coupling Vejj to a small (but nonzero) number; 
moreover, by small perturbations of v and/or q we can give Vf>ff any phase we 
please. But then, by putting a large number p of such chains in parallel, we 
can make the resulting v^jj lie anywhere in the complex plane we please. In 
particular, we can make Vejj equal to —q, which causes the partition function 
Zbs V to be zero. 



The proof of Theorem 2.2 is based on a complex- variables result due to 



Beraha, Kahane and Weiss |4C, 41, 42, 43 1; see [39] for details. 

3 A Bound in Terms of Maximum Degree 

Very little is known rigorously about the phase diagrams of Potts anti- 
ferromagnets, but one general result is available: for q large enough (how 
large depends on the lattice in question), the antiferromagnetic g-state Potts 
model has a unique infinite-volume Gibbs measure and exponential decay of 
correlations at all temperatures, including zero temperature. (Physically, the 
system is disordered as a result of the large ground-state entropy, so that 
zero temperature belongs to the high-temperature regime!) More precisely. 



using the Dobrushin uniqueness theorem [44, ^] it can be proven 1^^ that 
for a countable graph G in which every vertex has at most r nearest neigh- 
bors, the g-state Potts-model Gibbs measure on G is unique for all integer 
q > 2r whenever — 1 < fe ^ for all edges e. 

Uniqueness of the Gibbs measure means that first- order phase transitions 
are excluded, but higher-order phase transitions are not necessarily ruled 
out. In particular, the foregoing result does not imply the analyticity of the 
free energy for large-g Potts antiferromagnets, but it does make it plausible. 
Of course, a result that holds for integer q > qo need not hold for all real 
q > qo, much less for a complex neighborhood of that real semi-axis; but it 
does suggest that such a result might be true. This led me to conjecture that 
there might exist universal constants G{r) < oo such that, for all loopless 
graphs G of maximum degree < r, the zeros of the chromatic polynomial 
Pg{(i) [^-nd more generally of antiferromagnetic Potts partition functions] 
lie in the disc \q\ < C{r)!^ This is in fact the case, and the bound holds 
throughout the "complex antiferromagnetic regime" 1 1 + Ug | < 1 : 



Theorem 3.1 |14| There exist universal constants G{r) < 7.963907r such 
that, for all loopless graphs G = {V, E) of maximum degree < r, equipped 
with complex edge weights {ve}e£E satisfying |1 -|-?;e| < 1 for all e, the zeros 
of ZG{q,{ve}) all lie in the disc \q\ < G{r)vmax, where Vmax = maxlve]. In 

particular, the zeros of Pciq) all lie in the disc \q\ < C{r). 

This linear dependence on r is best possible, as the example of the complete 
graph Kr+i shows that G{r) > r. However, the constant 7.963907 can 
presumably be improved (see Section ^ below). 

It is amusing to note that the presence of one vertex of large degree 
cannot lead to large chromatic roots. More precisely, if all but one of the 
vertices of G have degree < r, then I can prove that the zeros of Pcio.) lie in 



the disc \q\ < C{r) + 1 |14|. Please note that a result of this kind cannot hold 
if "all but one" is replaced by "all but two", for in this case the chromatic 
roots can be unbounded, even when r = 2 and G is planar, as the graphs 
Qs,p show. 

The proofs of these results are based on well-known methods of mathe- 
matical statistical mechanics. The first step is to transform the Whitney- 
Tutte-Fortuin-Kasteleyn representation ( |1.4| ) into a gas of "polymers" inter- 
acting via a hard-core exclusion. One then invokes the Dobrushin condition 
1 49, ^ (or the closely related Kotecky-Preiss condition |51, ^) for the non- 



vanishing of a polymer-model partition function. Lastly, one verifies these 

^ I later learned that this result for Pg{<i) had been conjectured by mathematicians 
more than 25 years ago wA Mq, Question 6.1]. 
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conditions for our particular polymer model, using a series of simple combi- 
natorial lemmas. Details can be found in fl^ . 

With a little more work, it should be possible to extend the arguments 
of 1 14] to prove the existence and analyticity of the limiting free energy per 
unit volume (|1.5| ) for suitable regular lattices Goo and translation-invariant 
edge weights Ve, in the same region of complex q- and {ve}-space where Z 
is proven to be nonvanishing uniformly in the volume. In particular, this 
would provide a convergent expansion for the limiting free energy in powers 
of 1/q. However, I have not worked out the details. 

4 Some Conjectures and Open Questions 



The bound in Theorem |3.1| is, of course, far from sharp, and it is of 
some interest to speculate on what the best-possible results might be. Let 
us define 

Copt{f) = max{|g|: Pg{q) = for some loopless graph G of maximum degree r} 

(4.1) 
The example of the complete graph Kr+i shows that Coptir) > r. It is easy 
to see that Copi(l) = 1 and Copt{'2) = 2; and there is some evidence that 
Copt(3) = 3. But, at least for r > 4, Copt{r) must in fact be strictly larger 
than r, as is shown by numerical computations on the complete bipartite 
graph Kr^r [14|. Indeed, for large r it seems that Kr^r has a chromatic 



zero of magnitude approximately 1.5r. On the other hand, the Dobrushin 
uniqueness result quoted in the preceding section suggests that C{r) = 2r 
might suffice. 

We can pose these questions more generally as follows: Let Q he a class 
of finite graphs, and let V be a subset of the complex plane. Then we can 
ask about the sets 

Si{G,V) = U \J{q(^C: Zciq^v) =0} (4.2) 

Geg vev 

S2iG,V) = \J U {qGC: ZG(g,K}) = 0} (4.3) 

Gee {ve}: VeeVWe 

Among the interesting cases are the chromatic polynomials V = {—1}, the 
antiferromagnetic Potts models V = [—1,0], and the complex antiferromag- 
netic Potts models V = A = {v ^ C: \1 + v\ < 1}. Indeed, one moral of this 
work is that some questions concerning chromatic polynomials are most nat- 
urally studied in the more general context of antiferromagnetic or complex 
antiferromagnetic Potts models (with not-necessarily-equal edge weights). 
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The results of [0] show that the set S2{Gr,^) is bounded, where Gr is the 
set of all loopless graphs of maximum degree < r, and more generally that 
the set S2{Qr,A) is bounded, where Q'^. is the set of all loopless graphs of 
second-largest degree < r. But it would be interesting to examine in more 
detail the location of all these sets in the complex plane, and to prove sharper 
bounds. 

Another direction in which the results of [0] could be extended is by find- 
ing a criterion weaker than bounded maximum degree (or bounded second- 
largest degree) under which the zeros of Pg{q) and ZG{q,{ve}) could be 
shown to be bounded. An interesting idea was suggested very recently by 
Shrock and Tsai [H, |5^, who studied a variety of families of graphs and 
arrived at a conjecture that can be rephrased as follows: For G = {V, E) and 
X, y € y, define 

A(x,y) = max 7^ of edge-disjoint paths from x to y (4.4a) 

= min ^ of edges separating x from y (4.4b) 

and 

A(G) = maxA(2;,y) . (4.5) 

Clearly X{x,y) < min[deg(x),deg(y)] and hence A(G) < second-largest de- 
gree of G. Now let G^ be the set of all loopless graphs with A(G) < r. Then 
the conjecture is that the set S2{G^,V) is bounded, where V = {—1} or 
[—1,0] or perhaps even A. This possible connection of chromatic-polynomial 
and Potts-model problems with max-fiow problems is intriguing. Note that 
A(G) and A{G,{ve}) possess a "naturalness" property that maximum de- 
gree and its relatives lack: namely, for any graph G with blocks (2-connected 
components) Gi, . . . , Gb, we have A(G, {ve}) = max A(Gj, {ve})- 

l<i<b 
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